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We find a self-gravitating monopole and its black hole solution in Brans-Dicke (BD) theory. We 
mainly discuss the properties of these solutions in the Einstein frame and compare the solutions 
with those in general relativity (GR) on the following points. : From the field distributions of the 
generic type of self-gravitating monopole solutions, we can find that the YM potential and the Higgs 
field hardly depend on the BD parameter for most of the solution. There is the upper limit of the 
vacuum expectation value of the Higgs field to exist a solution as in GR. Since the BD scalar field 
has the effect of lessening an effective gauge charge, the upper limit in BD theory (in the lu — 
case) becomes about 30% larger than in GR. In some parameter ranges, there are two nontrivial 
solutions with the same mass, one of which can be regarded as the excited state of the other. This is 
confirmed by the analysis by catastrophe theory, which stated that the excited solution is unstable. 
We also find that the BD scalar field varies larger for solutions of smaller horizon radius, which can 
be understood from the differences of non-trivial structure outside the horizon. A scalar mass and 
thermodynamical properties of new solutions are also examined. Our analysis may give an insight 
into solutions in other theories of gravity; particularly, a theory with a dilaton field may show similar 
effects because of its coupling to a gauge field. 
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ON , After the discovery of a particlelike solution in the Einstein- Yang-Mills (EYM) system by Bartnik and Mckinnon 

jjj , a variety of particlelike solutions and black hole solutions with non- Abelian fields were found B-|lql . Some of 

them are the solutions in which the self-gravity is taken into account for the solution in flat space-time. The others 

O" 1 are new types of solutions which can exist in the presence of gravity. They show interesting properties, which can not 

; I \ be seen in the flat spacetime solutions or well known Kerr-Newman black hole, such as stability, thermodynamical 

bJQ' properties, mass inflation phenomena and so on |14j. In the context of a counterexample for the black hole no-hair 

conjecture, a monopole black hole PT13J may be the most important among them, since both the Reissner-Nortstrom 

• ^h , (RN) black hole and the monopole black hole with the same mass and the gauge charge are stable within a certain 

parameter region and it is constructed in the EYM-Higgs (EYMH) system which is not a phenomenological model 

such as the Skyrme and the Proca model. The monopole black hole solution describes a space-time structure with 



a 



I. INTRODUCTION 



an event horizon in the gravitating 't Hooft-Polyakov magnetic monopole [la]. The existence of such a black hole 
solution was first expected by examining a simple distribution of the energy density, where it is almost constant inside 
a core and decays as r~ 4 outside ||. Later Breitcnlohner et al. JO]], Aichelburg |12| and we [|13| found the "exact" 
solutions numerically. 

In the big bang scenario, the 't Hooft-Polyakov monopole is overproduced during phase transitions in the early 
universe. This is one of the reasons why the inflationary scenario was proposed Jig ]. Although such an inflationary 
scenario was discussed originally in general relativity (GR), it turns out that introduction of a scalar field (such as 
the Brans-Dicke (BD) scalar field) can make a great change in scenarios of the very early universe U%. Moreover, 
Vilenkin and Linde proposed that inflation occurs inside the 't Hooft-Polyakov monopole |lq] , where no fine tuning 
of the initial value is necessary. Such a monopole inflation is also changed qualitatively by introducing a scalar field 
Jl9| . So properties of the gravitating monopole and its black hole solution may be greatly modified in scalar-tensor 
theory compared with in GR. 
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As for the dilatonic black holes in an effective theory of superstring, it is known that such a scalar field also affects 
many features [J20|]. In BD theory, it is also known that the Kerr or Kerr-Newman black hole is a unique solution 
in the vacuum case or in the case with the Maxwell field, because of the conformal invariance of the Maxwell field 
and the black hole uniqueness theorem f2l]| . Recently, however, we found non-trivial black hole solutions with the 
non-Abclian field in BD theory E2] . They are generalizations of neutral non- Abelian black holes in GR (Proca ||] and 
Skyrme black holes || ) . Then, it is important to study a self-gravitating monopole and its black hole solution in BD 
theory as the charged solution. In addition, we have the following reasons to consider the monopole black hole in BD 
theory, (i) For the neutral non- Abelian black holes, we can not find a more stable solution than a trivial solution, i.e., 
Schwarzschild solution even in BD theory. We expect that the monopole black hole can be more stable than a trivial 
solution, i.e., Reissner-Nortstrom (RN) solution in the BD-YMH system. So it would be more important than the 
neutral type. Particularly, when we discuss the Hawking radiation pjl of the RN black hole in the EYMH system, it 
may evolve into the monopole black hole and eventually into the gravitating monopole, which would be a candidate 
for its remnant. The scalar field may cause a serious effect on this evolution, (ii) To investigate the stability of the 
neutral non- Abelian black holes, we can apply catastrophe theory even in BD theory ||J2^]. Using the variables in the 
Einstein frame, we find that the BD scalar field is less important from the perspective of catastrophe theory. However 
the solution is not evident in the charged monopole black hole case. 

This paper is organized as follows. We introduce basic Ansatze and the field equations in the BD-YMH system in 
Sec. II. In Sec. Ill, we present numerical solutions for the regular self-gravitating monopole. In order to see their 
structure, we show their field configurations and compare them with the solution in GR. We also discuss the upper 
limit of the vacuum expectation value of the Higgs field for a solution to exist. In Sec. IV, we present the monopole 
black hole solution. In some parameter regions, there exist two types of non-trivial solutions, one of which is stable 
and the other is not. We show that the dependence of field configurations on the BD parameter is different in each 
type. This shows the difference of the black hole structures which can not be seen only by considering them in GR. 
We also discuss upper limits for the vacuum expectation value of the Higgs field for a static solution to exist in black 
hole cases. The concluding remarks will follow in Sec. V. Throughout this paper we use units of c=H=l. Notations 
and definitions follow Misner-Thorne- Wheeler fl24|. 



II. BASIC EQUATIONS 



The action of BD theory is written as J2E 



JL(^R_^V <>Y', 



L r . 



(2.1) 



where k 2 = 8irG with G being Newton's gravitational constant. u> is the BD parameter, and L m is Lagrangian of the 
matter fields. The dimensionless BD scalar field <f> is normalized by G. In our previous paper p2[ , we investigated 
two types of non- Abelian black hole solutions in BD theory (Proca black hole and Skyrme black hole) and found 
that there are several advantages to discussion in the Einstein conformal frame rather than in the BD frame. For 
example, the definition of thermodynamical variables of the solutions and the stability analysis using catastrophe 
theory. Although, of course, we can reach the same results in the BD frame, we transform the Einstein frame Cfl] as 
a matter of convenience. The conformal transformation in the present model is written by 
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Then, we find the equivalent action S = S/(f>o given as [t 
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where 
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Variables with a caret " denote those in the Einstein frame. L m is different from L m unless the matter field is 
conformally invariant, as with the Yang-Mills field. The constant 0o = 2(2 + cu) / (2oj + 3) is the value of the BD scalar 
field <j) a t spatial infinity, which guarantees that we can regard G as Newton's gravitational constant |25|. In this 
paper we investigate the YMH system, where Lagrangian L m is written as 
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F a is the field strength of the SU(2) YM field and expressed by its potential A a as 

F% = dpA^ - d v A\ + et ahc A\A\, 
with the gauge coupling constant e. $ a is the real triplet Higgs field and D M is the covariant derivative: 

D^ a = d^ a + ee abc A b ^ c . 



(2.6) 
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The theoretical parameters v and A are the vacuum expectation value and the self-coupling constant of the Higgs 
field, respectively. Below we consider them only in the Einstein conformal frame, and drop the caret ". 
We assume that a space-time is static and spherically symmetric, in which case the metric is written as 
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For the matter fields, we adopt the hedgehog ansatz given by 

$ Q = vr a h(r), 



(2.10) 



A a = 0, 



Ap=w li e r . (^ = 1,2,3), 



where r a and w c „ are a un it ra dial vector in the internal spa ce a nd a triad, respectively. 
Variation of the action (2.3) with the matter Lagrangian (|2.6|) leads to the field equations 
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r ~ r _(l-w 2 Y , v -to(uh\ , A _ 2 ^» a 
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We have introduced the following dimensionless variables: 

f = et>r, m = Gevm, <p = ncp, (2-21) 

and dimensionless parameters: 

i = vVG, A = A/e 2 . (2.22) 

Although the solution exists when v ~ Mpi, where Mpi is Planck mass, it can be described by a classical field 

configuration in the limit of a weak gauge coupling constant e, because its Compton wavelength ~ e/v is much 

smaller than the radius of the classical monopole solution ~ 1/ev in this case. Moreover, since the energy density 
~ e 2 v 4 <§; M p[ , we can treat the system classically. 
The boundary conditions at spatial infinity are 

lim m = M < oo, (2.23) 

lim 6 = 0, (2.24) 
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lim k/3lp = lim In -f- = 0, (2.25) 



lim h=l, (2.26) 
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lim w = 0. (2.27) 
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These conditions imply that the space-time approaches a flat one and that the solution is a charged object. The 
field equations are singular at the origin and at the event horizon for a self-gravitating monopole and its black hole 
solution, respectively. 

For a self-gravitating monopole solution, we require regularity at the origin. Expanding the Eqs. ( |2.13 )-(2.17) 
around the origin, we find that regularity is guaranteed when the fields behaves as 

„,( e ) « l + IcU 2 + 0(e 3 ), (2.28) 

h(e)^C h e + 0(e 3 ), (2.29) 



<p(e) « C v - ^e-^v 2 Cl + ^e"^ e 2 + 0(e 3 ), (2.30) 



5{ e)~C & -^(cl + e ^Ct)r- ()u :i ). (2.31) 



m(e) « -^ ( C£ + e-^Cl + V 2 ^ ] 6 3 + 0(e 4 ). (2.32) 

0o V 6 / 



C w , Ch and C v are constant, which should be determined iteratively so that the boundary conditions (2.25)-( 2.27 ) 
are satisfied. We also require that no singularity exists, i.e., 
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In the case of the monopole black hole, we assume the existence of a regular event horizon at r = Th- For the 
metric 



(2.34) 
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For the matter fields, the square brackets in Eqs. ( 2.15 )-( 2.17 ) must vanish at the horizon. Hence we find that 
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We also require that no singularity exists outside the horizon, i.e., 
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Hence, we should determine the values of C w , Ch and C v (self-gravitating monopole case) or ipu, Whi and ha 
(monopole black hole case) iteratively so that the boundary conditions at infinity are fulfilled. In this sense these are 
shooting parameters. It is a remarkable property of the nontrivial structure of many systems discussed before that we 
can not choose arbitrarily the values both of the fields and of their derivatives at the horizon. When we investigate 
the internal structure of non-Abelian black holes, we may expect the e xi stenc e of the inner (Cauchy) horizon. There, 
the fields must satisfy the same type of constraint such as Eqs. ( 2.36 )-( pT38 ). However, the shooting parameters on 
the black hole horizon which satisfy the boundary condition at infinity do not necessarily fulfill the constraint at the 
inner horizon. By this behavior the mass function diverges and a mass inflation phenomenon is realized inside of the 
black hole event horizon for non-Abelian black holes vM - 

A non-trivial solution does not necessarily exist for the given theoretical parameters. However, for arbitrary values 
of v and A, there exists a RN black hole solution such as 
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M is the gravitational mass at spatial infinity and Q = 2v^/tt/(J)o is the magnetic charge of the black hole. The radius 
of the event horizon of the RN black hole is constrained to be fn _ Q- The equal sign case implies an extreme solution. 
If we equate fi = 0, (i.e., u> — > oo) the model recovers the EYMH system and we find a self-gravitating monopole and 
its black hole solution investigated in Ref. p-fi3fl. In addition, if we put G = 0, i.e., neglect the gravitational effect, 
the model recovers the YMH system, which has a famous 't Hooft-Polyakov solution for a regular magnetic monopole 
[Ell. Further, assuming A = 0, we obtain the BPS monopole solution pq|, 
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w(f) = —— — -, h(f) = cothf — — . 
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In the next section, we show the solutions of a self-gravitating monopole both in GR and BD theory. 
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III. SELF-GRAVITATING MONOPOLE IN BRANS-DICKE THEORY 



The self-gravitating monopole in GR were investigated in detail by several authors |9|-|13| and found to have many 
interesting properties, some of which can not be seen for the monopole solution in flat space-time. We list their 
properties. (1) The distribution of the non-Abelian structure decays exponentially with respect to the distance from 
the origin and the monopole has its core radius ~ 1/ev. (2) The solutions are characterized by the node number of the 
YM potential, which is different the from topological winding number of the monopole. (3) There exists the maximum 
parameter v max beyond which there is no monopole solution. v max weakly depends on A. (4) There exists the critical 
parameter v extr eme for the small A. At v — v ex treme, the monopole solution changes to extreme RN solution. Between 
v extreme and v maXl two monopole solutions with a different mass appear for each value of v. The higher mass solution 
is unstable while the lower one is stable. v ex treme also weakly depends on A. 

Now we turn to the monopole solutions in BD theory and compare their properties with those in GR from the above 



points of view. Solving the system (2.13)-(2.20) with the boundary conditions (2.23)-(2.33), we obtain gravitating 
monopole solutions. Although they are found only numerically, there exists a rigorous proof of the existence of self- 
gravitating monopole and its black hole solution in the A — » oo case in GR |11[. We expect that the same discussion 
holds in the more general case including BD theory. 

We show the field configurations of the gravitating monopole in Figs. 1. We chose the parameters as v = O.f , 
A = 0.1 and several values of u) = —1.4, —1.0, 0, oo. The oj = oo case corresponds to GR. (a) is r-h,w diagram. We 
can see that the configurations of the YM potential and the Higgs field hardly depend on the BD parameter u>. This 
implies that the solution can have almost the same structure as the 't Hooft-Polyakov monopole in flat space-time, 
where the repulsive force of the gauge field balances with the attractive force of the Higgs field, (b) and (c) are r-m g 
and r-S diagram, respectively. Since fh g is defined as 
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it diverges at the origin. The mass function rh g strongly depends on ui. For a different BD parameter each solution 
has a different effective global gauge charge Q. For the small u>, the charge is small, i.e., the monopole in BD theory 
has a smaller charge than that in GR. Since there is a prefactor Q in the S and m equations, it contributes directly 
to the behavior of the metric functions. There is another factor which may change the behavior of metric functions. 
It is the BD scalar field p. However from Fig. 1 (d), <^3<0.07 even for u> — —1.4, which implies the e _/3 ^ rs 1. dp/df 
term is also smaller than other terms in the field equations. Hence we can neglect the dp contribution to other fields 
in the zeroth approximation. However, we may not neglect the dp contribution for the parameter range v rs v ex treme 
discussed later, where the gravitational effect plays an important roll. 

In Fig. 2 we show the distribution of the energy density p = — Tq for the same parameters as Fig. 1 (It is normalized 
as p = pG/(ev) 2 ). An arrow shows the core radius of the monopole. The prefactor l/4>o for the matter Lagrangian 
can be absorbed by defining the new variables and the effective parameters as A a = A't/yfcpQ, $ = $/y^o, e = \/0o e > 
v = v/y/ifo and A = 4>q\. Hence the core radius ~ 1/ev = 1/ev hardly depends on (f>o, i.e. BD parameter lu. The 
energy density is almost constant inside the monopole core, but it decays as 1/f 4 outside. Note that although w 
decays exponentially, the YM potential decays ~ 1/r for a large radius. 

As for the monopole in flat space-time, the uniqueness of the solution was proven p9| , which means there is no 
excited mode of this solution. As for gravitating counterpart, however, there exist radial excited solutions characterized 
by the node number n of the YM potential for small v in GR [O . We also find such excited solutions in BD theory. 
In the v — > limit, the Higgs field becomes constant and the YM field recovers its apparent gauge symmetry. Then 
the YM field is conformally invariant and the BD scalar field becomes trivial by the regularity. As a result, the n = 1 
excited solution approaches a regular Bartnik-McKinnon solution [nj in the EYM system in this limit as discussed in 
Rcf. [[llj. Hereafter, we focus on the solution with n = 0, i.e., the ground state solution. 

In GR, there is the maximum value of the vacuum expectation value v max for a fixed A where a regular monopole 
solution can exist O]. We can understand this intuitively as follows. We compare the mass of the monopole ~ Awv/e 
with its core radius ~ 1/ev. When v becomes large, the mass of the monopole increases, as does the gravitational 
radius ~ 8nGv/e, while its core radius shrinks. Hence, its core radius eventually gets into gravitational radius and the 
solution may develop into the RN black hole being swallowed in its non-trivial structure by the horizon. The existence 
of the boundary of the parameter is interesting also in the context of the inflation. Since there is no non-trivial static 
configuration beyond the maximum parameter, the solution must take a configuration of the static RN black hole or 
a nonstatic configuration. One of the possibilities in the latter case is the inflating solution. Actually it is confirmed 
that the boundary of the static solution almost coincides with the parameter region where the monopole inflation 
occurs [1811. 



In BD theory, it was shown that when we consider the global monopole, any inflating monopole eventually shrinks 
contrary to the case in GR [|19|. Although we consider gauge monopole here, the BD scalar field may give a serious 
effect. Moreover it is found that the gravitating monopole solution in TV = 4 supersymmetric low-energy superstring 
theory exists for arbitrary values of vacuum expectation value of the Higgs field, where the dilatonic field plays an 
important roll [p0[. Sin ce t he BD scalar field can be regarded as the dilaton field in a special coupling constant case 
as can be seen in Eq. ( |2.3| ), there is a possibility that v max disappears. Actually, taking v larger, we find the BD 
scalar field gives a larger contribution. Hence we are interested in examining the maximum value v max in BD theory. 

First, we discuss analytically the parameter region where the nontrivial solutions exist following the analysis in 
Ref. M. For the monopole solution, the mass function behaves as Eq. ( [2.32 ) around the origin by regularity. At 
spatial infinity, the solution approaches RN solution, 
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fh(f) =M- — - + • • • . (3.2) 
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From this boundary condition, g rr should have a minimum which corresponds to the maximum of the m/r. This 
occurs around r m ~ ^irv 2 / foM . For the monopole in flat space-time, we have 

Attv 2 
Mfiat = -—5(A) , (3.3) 

fo 

g(X) takes values from 1 to 1.787 when A changes from to oo. We assume M oc v 2 even for the gravitating monopole 
in both theories. Under this assumption, r m ~ 1/g ~ 1. So g rr takes a minimum value of about g rr ~ 1 — 4irv 2 / fo. 
When v — v extreme = \/fo/4ir, g rr takes and the degenerate horizon appears at r = r extreme . Because the region 
f < f extreme is separated infinitely in geodesic distance from the outer region |31| , the monopole structure could be 
thought to be confined in the region r < f extreme and the solution looks like a RN black hole with the horizon radius 
r e xtreme- Naively Vextreme gives the maximum value of v. As a result, we can expect that the critical value v ex treme 
becomes large for the large fo, i.e., the small value of to. 

Breitenlohner et al. confirmed the above discussion numerically in GR (fo — 1) and found the maximum value v m ax 
for the fixed A by examining the minimum of g rr [0| . Moreover they found the parameter region of A where v 'extreme is 
not the maximum value of v. In the parameter region v extreme < v < v max , there are two different nontrivial solutions 
with the same v which is not obvious in the above simple analysis. On v-M diagram, these solutions construct two 
solution branches divided by a cusp structure at v = v max . The lower mass branch ends up with v max and the higher 
mass branch is connected to the RN black hole branch at v = v extreme- The cusp structure is important for the 
stability analysis by using catastrophe theory |13j. For the larger A, v max and v cr it coincide [p3| . Since v max takes a 
value very close to v extreme, we search v extreme by the behavior of g rr first and estimate v max . 

In the same manner, we examine v max in BD theory with lu — 0. Figs. 3 shows the field configurations of the 
monopole with A = 0.2 and around the maximum parameter v max ((a) f-g rr (b) r-h, w). We can find that g rr has 
the minimum around r = r extreme ~ 1-53 and the minimum value decreases as v becomes larger. When v approaches 
Vmax , the behavior of the YM potential and the Higgs field changes more rapidly near the f extreme and takes w ~ 
and kfsl, i.e., the solution can be approximated by a RN black hole solution outside of f extreme- When we have 
Vextreme, the horizon degenerates and the outer region becomes the extreme RN black hole. Since the gravitational 
effect becomes strong around this parameter region as seen from Fig. 3 (a), the amplitude of tp and the dependence 
of w and h on the BD parameter can not be negligible. We show the v m ax-^ dia gram w ith 10 — and oo in Fig. 4 
(the solid lines). Though v max is larger compared to the estimation Vextreme = \ffo[^ ~ 0.282 (for GR), ~ 0.325 
(for u = 0) due to the naive estimation, we can obtain the relevant dependence of v max on fo (or to). If we take a 
more accurate value of r m calculated numerically, we have the better estimation of Vextreme ~ 0.369 (for GR) and 
Vextreme ~ 0.497 (for lu = 0) As we describe later, the similar phenomena can be seen in monopole black hole. These 
are also shown in this diagram by the dotted lines (In this diagram, f# = O.4.). 

IV. MONOPOLE BLACK HOLE IN BRANS-DICKE THEORY 

The monopole black hole solution in GR has the following properties in addition to (l)-(2) of Sec. III. The number 
of the solutions depends on the parameters complicatedly, since another parameter, the horizon radius, is added in 
the black hole case. We show the schematic diagram of v-ru for small A in Fig. 5 [B9[. In the shadowed region, the 
nontrivial solution exists. The f# = axis corresponds to the regular monopole solution. (5) For a fixed value of 
Th < ?*h (See Fig. 5), we find two characteristic values v max and v ex treme- They are extensions of properties (3) and 
(4) in the regular monopole case discussed in the previous section. In the region (i) in Fig. 5 (v ex treme < v < V max ), 



there are two nontrivial solutions with different mass. The solution branch with higher mass is connected to the 
extreme RN black hole solution in the v-M diagram although the horizon radius changes discontinuously. (6) The 
solution in the lower mass branch is stable while the higher one is unstable. (7) For a fixed value of v < v* (See 
Fig. 5), there are two characteristic values rH,max an d rn^RN- In the region (ii) in Fig. 5 (rH,RN < r < rn.max), 
there are two nontrivial solutions. The solution branch with lower entropy (smaller fjj) is connected to the RN black 
hole at fjj = fn,RN- (8) The solution in the high entropy branch is stable while the lower one is unstable. (9) The 

regions (i) and (ii) disappear for A > A^L ~ 0.287 and A > \™ it ~ 0.696, respectively. (10) When we consider the 
thermodynamics of the monopole black hole, its temperature diverges in the minimum mass limit unlike the RN black 
hole. Through a evaporation process from a large RN black hole, the solution experiences the first and/or second 
order phase transitions. 

We turn to th e mono pole black h ole in BD theory. Solving the system ( 2.13 )-( 2.20| ) numerically with the boundary 



conditions (2.23)-(2.27) and (2.34)-(2.40), we obtain a black hole solution. Figs. 6 are the field configurations with the 
parameters A = 0.1, v — 0.1, u) — and the horizon radius ?h — 0.2, 0.4 and 0.6 ((a) r-h,w (b) f-rh g (c) f-5 (d) f-<p) 
We also plot the gravitating monopole solution, i.e., fa = 0, for comparison. As the horizon radius becomes large, 
the curvature of the space-time outside the horizon would become small resulting in the variation of the BD scalar 
field being small. Hence the difference between the solutions in BD theory and in GR appears most conspicuously 
in the monopole solution case. The boundary value of w and h approaches w(rn) — > and h(rn) — > 1 for a large 
horizon radius. Finally the solution becomes a RN black hole in the rn = ?h,rn limit. The behaviors of the solutions 
around this limit depend on the parameters complicatedly as discussed later. 

In Fig. 7(a), we show the M-th diagram for A = 0.1, v — 0.1 and oj = —1.4, — 1, and oo. Dotted lines show RN 
solution branches. Note that the RN solution in BD theory is different from that in GR because it has a different 
effective gauge charge by the factor cf) . Unlike the RN black hole, the monopole black hole does not have a extreme 
limit but has the fjj = limit, which corresponds to the gravitating monopole solution. Fig. 7(b) is a magnification 
around the maximum horizon radius of each theory. As is the case of GR, there are two solutions in some ranges of 
horizon radius. This parameter region corresponds to the region (ii) in Fig. 5. We call the upper (lower) branch a 
high (low) entropy branch. In both theories, these branches form a cusp structure at the point A, where the solution 
has the maximum horizon radius rH,max- At the point B, the solution coincides with the RN black hole solution with 
the horizon radius th.rn- The solutions have a qualitatively similar structure in both theories. The maximum mass 
and the maximum horizon radius decrease when u> becomes small. Furthermore we find that the range of horizon 
radius where the low entropy branch exists becomes larger when we take the w as small. These behaviors are mainly 
due to the effect of the boundary value of 4>q. As we mentioned, v — v/\Q>o can be considered as an effective vacuum 
expectation value in BD theory. Since it becomes small when u) becomes small, we find that the width of region (ii) 
in Fig. 5 becomes large. 

Next we consider a difference of solution in each branch. We consider the solution with v = 0.1, A = 0.1 and 
f/j = 0.59, and change the BD parameter to. Figs. 8 and 9 show configurations of the fields ((a) f-h, w (b) f-fh g (c) 
f-(p ). In the high entropy branch (Fig. 9), the YM potential and the Higgs field hardly depend on uo. This is the same 
as the self-gravitating monopole case and the structure is supported by the balance between these fields. On the other 
hand, they have small amplitude in the low entropy branch (Fig. 10). The Higgs field can be approximated as h « 1, 
which means that it takes almost its vacuum value, while w « 0, where the YM potential decays as f . Hence, what 
attracts the YM field is not the Higgs field but the self-gravitational attractive force. Actually, the counterpart of 
the low entropy branch does not exist in flat space-time. The dependence on lo in the low entropy branch is due to 
the difference of the effective gauge charge again. The horizon radius fh = 0.59 line crosses the low entropy branch 
relatively near the point A (See Fig. 8(b)). In the large u> (GR), however, it crosses around the point B, where the 
solution coincides with the RN black hole. Hence the field configurations approach a trivial solution. 

In the present model, the YM field obtains the mass through a spontaneous symmetry breaking mechanism and 
the trace part of the energy-momentum tenser does not vanish. In this case the BD scalar field takes a nontrivial 
configuration and we can define the scalar mass M s by the asymptotic behavior of the BD scalar field as 

l + ™±). (4.1) 



Fig. 10 shows uj dependence of the gravitational mass M and scalar mass M s for fixed horizon radius f/j = 0.59. As 
the solution can be approximated by the RN black hole in the zeroth approximation, M behaves as 



, ~i 2ttw 2 
Attv 



2f h 
This gives good agreement with the u dependence in Fig. 10(a). M s can be expressed as 



(4.2) 



*• = 5Jrl*> (4 - 3) 

where s is the sensitivity of a self-gravitating object. Though s varies as ix> changes, its configuration is negligible for 



large w and M s behaves as ~ 1/w. Since s — » 1/2 in the w — > —3/2 limit, the fraction in Eq. (4.2) is indeterminate. 
However, Fig. 11(b) shows that M s takes a non-zero finite value. As M takes almost the same value in each branch, 
the solution in the high entropy branch has smaller sensitivity than the lower branch. 

In GR, it is considered that when A becomes large, the parameter point where the Fh.rn = ^H.max in Fig. 5s shifts 
to the left side, i.e., smaller v and the region (ii) becomes small. When A takes the critical value \™lf ^ ne re gi° n (ii) 
finally disappears |ll],|4| . Although we do not examine the definite value of a^u m BD theory, the dependence on uj 
would be small because A = A/e 2 = A/e 2 does not depend on cj) . 

The stability of these exotic objects is one of the main interests. Aichelburg [O] analyzed the monopole black hole 
and its stability via the linear perturbation method in GR and showed that most of the solutions are stable within 
the limit of A — > oo. They also obtained the solutions in the finite A case and suggested that if we have two monopole 
black hole solutions for a fixed value of v, the low energy one is stable while the other is not. Similar works arc found 
in |llj, where the authors considered the weak coupling limit. We examined stabilities via catastrophe theory |L3| and 
obtained a unified picture of the stabilities of the gravitating monopole and the monopole black hole. In catastrophe 
theory [|40j, we discuss the shape of a potential function as the variation of the control parameters. In that paper, we 
adopted A, v and M as the control parameters, entropy S as the potential function and 5{rh) as the state variable, 
and found that the system is classified into a swallow tail catastrophe. In general, a swallow tail catastrophe has three 
control parameters a, b, c and one state variable x and its potential function is described as 



We regard the parameters as 



V^ ' ic r*V\f o in^iH Ktt- n\ \' ' 



V = 7 x 5 + -ax 3 + -bx 2 + ex . (4.4) 

Do Zi 



a = a(\,v), b = b(M,v), c = 0. (4.5) 



KirL 1S obtained by a(X cr l t ,0) = 0. Although we have a new control parameter m in BD theory, the situation does 



not change. By analyzing the dependence of the solution on the parameters, Eq. (4.5) is changed to 



a = a(\, v,uj), b = b(M, v,lj), c = 0, (4-6) 

i.e., the BD parameter is not an intrinsic parameter but a dummy parameter in catastrophe theory, and does not 
result in qualitative difference from GR. This is the same as in the neutral type non- Abelian black hole case p2| . As 
a result, the monopole black hole in BD theory is also classified as a swallow tail catastrophe. From this we find that 
the solution in the high entropy branch and RN black hole with larger horizon radius than th,rn are stable, while 
the solutions in the low entropy branch and RN black hole with smaller horizon radius than rH,RN are unstable. 

These properties are interesting in the context of black hole no-hair conjecture and the uniqueness of the solution. 
For A > Xcrlv t nere are two solutions (monopole black hole and RN black hole) when we fix the global charge M 
and Q. This is the counterexample of the strong no-hair conjecture, i.e., if the global charges are fixed, there exists 
only one solution. Moreover when A is smaller than A^? it , there is a mass range where two stable solutions (monopole 
black hole in the high entropy branch and RN black hole) exist for a fixed mass. This violates even the weak no-hair 
conjecture which states that if the global charges are fixed, there exists only one stable solution. These behaviors are 
shared by both theories. In GR, however, we have no idea how to identify whether a black hole is a monopole black 
hole or RN black hole from infinity, because the NA structure and the Higgs field decrease exponentially. On the 
other hand, in BD theory we can obtain the information about the scalar mass, which can be used to tell the type of 
black holes. Hence we can identify each black hole uniquely from infinity. 

When we take into account the quantum effect, a black hole loses its energy by the Hawking evaporation process. 
The evolution of the monopole black hole by this process is rather complicated. First, we show an inverse temperature 
in the A = 1.0, v = 0.1 and to — 0, oo (GR) in Fig. 11(a). The inverse temperature of the RN black hole has the 
minimum \jT = 6^/3ttQ at M = 2Q/V3 and diverges at the extreme limit M = Q, where the evaporation ends. On 
the other hand, the inverse temperature of the monopole black hole becomes zero in the M — M monopo i e limit as with 
the Schwarzschild black hole. Consider the evolution of a RN black hole with fu > ?H,max by evaporation process. In 
GR in this figure, it loses mass, and when the solution approaches the minimum of the inverse temperature, its heat 
capacity changes discontinuously. [Note that its specific heat changes its sign continuously] This is the second order 
phase transition. The black hole loses its mass further and reaches the point B. Below this mass, a RN black hole 



becomes unstable and the solution traces the monopole black hole branch. The nontrivial structure comes out of the 
event horizon gradually. This is also the second phase transition. At this point, however, both the heat capacity and 
the specific heat change discontinuously. After that, the evaporation continues until the event horizon disappears. 
Finally, the regular monopole may remain as a remnant. For u) = 0, however, it experiences a different evolution, 
because before it reaches the minimum of the inverse temperature, it reaches the point B and traces the monopole 
black hole branch continuously. 

Furthermore, the A < A^rit case ^ s different. Fig. 11(b) shows the inverse temperature of the solutions with the 
A = 0.1, v = 0.1. In the ui = —1.4, —1.0 and case, the RN black hole evolves to the point B as in the above case. 
However, since the RN black hole becomes unstable, the solution jumps to the stable monopole branch, increasing 
its entropy and temperature discontinuously (arrows in Fig. 11(b)). This is the first order phase transition. The 
nontrivial structure comes out suddenly (or catastrophecally) from the horizon. In the u> — oo case, the solution 
experiences the second order phase transition at M — 2Q/^/S before the first order phase transition. As shown above, 
the type and number of the phase transitions depend on the physical parameters. This property is more important 
for the stability of the solution in the system surrounded by a heat bath. In that system, stability change occurs at 
the point where the heat capacity diverges |4l[ . Another thing that we should note is that the monopole black hole 
objects are forced to be produced by this process if the accretion to the black hole is very little. 

Until now, we have studied the solution in the rather small value of v corresponding to the left side of Fig. 5. and 
found the region (ii) in BD theory. Next we consider the v max in BD theory. By the same method in the gravitating 
monopole case, we found v max . Fig. 4 shows the existence region of the black hole solution (fa — 0.4) in v-X plane 
for lu — and oo. In the left region of each dotted line, there are black hole solutions. When we take A as larger, 
the dependence of v m ax becomes smaller. We can not see the intrinsic difference also in this diagram. It is different 
to clarify the detailed structure of region (i) in BD theory because of the fine tuning of the asymptotic value of the 

BD field. In GR the region (i) disappears at A = X^ it ~ 0.287. As discussed before, since A does not depend on </> , 
A") would not be changed much even in BD theory. However, since the gravitational effect becomes strong around 
v = v max , the BD scalar field may affect the value of X cl : it . It needs further analysis to say definite thing. 

V. CONCLUSION 

We found the self-gravitating monopole and its black hole solution of monopole in BD theory and compared them 
with those in GR. We clarified the following important aspects. The configuration of the YM potential and the Higgs 
fields of a self-gravitating monopole and a monopole black hole hardly depend on the BD parameter in most of the 
mass range, which shows the structure resembles the 't Hooft-Polyakov monopole. We found the parameter region 
(ii) in Fig. 5, where two monopole black hole solutions exist for fixed horizon radius. Each solution branch has a 
different dependence for the BD parameter. For the higher entropy branch, the configurations of the fields hardly 
depend on the BD parameter, while the configuration of the fields changes to those of RN black hole in the v — > v max 
limit for the lower entropy branch. However, we can not find the region (i) because fine tuning is needed to satisfy 
the boundary condition of the BD field at infinity. One of the characteristic properties in BD theory is lessening 
the effective global gauge charge. By this effect, the M-fh diagram (Fig. 7) shifts to the left side, i.e., the mass 
of the solution becomes small when u> gets small. From this diagram, we can derive information about stability by 
applying catastrophe theory. As in the GR case, the high (low) entropy branch is stable (unstable). And we found 
that the BD parameter is not an intrinsic parameter but a dummy parameter in catastrophe theory. The effective 
global charge also affects the configuration for the high entropy branch, and does affect the existence region of the 
solutions in parameter regions. We investigated a boundary value of v above which there exist no static solution for 
fh = 0.4 and a self-gravitating solution. It is interesting that the allowed parameter region is extended in BD theory. 
We also discuss the thermodynamical properties of the black hole solutions. The number and the order of the phase 
transitions depend complicatedly on the theoretical parameters. 

These aspects may give us some important information about monopole black holes in other theories of gravity. In 
particular, we can anticipate the results in dilaton gravity because of its similarity to BD theory. Though we can not 
see the intrinsic qualitative difference form GR, some effective theory of unified theory may cause important effects 
for monopole fl30[ ] or its black hole solution, which may predict some information to the remnant of the Hawking 
evaporation. 

There are some issues which we should comment on. In this paper, we transformed the system from the original 
BD frame to the Einstein frame and discussed only in the latter frame. However, the solutions in the former frame 



are easily obtained from our solution by the inverse transformation of Eq. (2.2). For example, M-th diagram in the 



BD frame is obtained as follows. In the BD frame, a gravitational mass Mbd is defined by use of the scalar mass as 
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Mbd = M + M s . (5.1) 

The Mbd-Th diagram is shown in Fig. 12, where we normalize as Mbd = GevMbd- Because of the existence of a 
scalar field, a cusp structure disappears at the points A. This may show that Mbd is not appropriate to a control 
parameter of catastrophe theory [g2[Q. Since the solution in a high entropy branch has a larger scalar mass than 
that in a low entropy branch as shown in Fig. 10, the high entropy branch slides to the right of the lower branch in 
the BD frame. 

We note the relation between our results and other topological defects such as a global monopole or a cosmic string 
in BD theory or in more general scalar-tensor theories. For global monopole, it is investigated that a massless dilaton 
field or the BD scalar field affects its structure, and is discussed that there is no solution except the special value of 
the coupling constant p5fl . Similar arguments are seen for a cosmic string in ]3q ] where they discussed that "wire 
approximation" which can be valid in GR will not be used. For a massive dilaton field, though both global monopole 
and cosmic strings would not be affected seriously compared to the massless case, it is discussed that it enables us 
to decide the dilaton mass from observation |37|| . As for the texture, though exact results are almost unknown but 
Dando using self-similar ansatz showed that regular solution exists even in a massless dilaton case for a suitable choice 
of coupling constant |J38fl . Contrary to these cases, both the self-gravitating monopole and its black hole solution exist 
for any value of the BD parameter in our model. This would be due to the asymptotic mild behavior of the BD scalar 
field. 

We finally note that much work on the monopole solution has been developing from the context of superstring 
theory [|45{ . Recently monopole solutions in flat space-time, self-gravitating monopole solutions and its black hole 
solutions were obtained in Born-Infeld type action [E6j . However, little has been clarified about them. It should be 
interesting to investigate their structure, especially the internal structure around the central singularity, because of 
the remarkable form of its action. 
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FIG. 1. Field configurations of the self-gravitating monopole for v — 0.1 and A = 0.1 ((a) f-h,w (b) f-m g (c) r-8 (d) f-(p). 
We compared these distributions for u> = —1.4, — 1, 0, oo. Though metric functions and the scalar field depend on the BD 
parameter mainly due to the difference of the effective global charge, we can find that YM field and the Higgs field hardly do. 



la. 




FIG. 2. The density distribution of the monopole solution for the same parameters as Fig. 1. An arrow shows the core 
radius of the monopole which hardly depends on the BD parameter. The energy density decays as 1/f 4 outside the core. 
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FIG. 3. Field configurations of the self-gravitating monopole with A = 0.2 and uo — near the v m ax ((a) f-g Tr (b) r-h, w). 
Outside the minimum of g rr , the solution can be approximated by RN black hole. 
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FIG. 4. The allowed region of the static solution for the self-gravitating monopole (the solid lines) and its black hole with 
fn = 0.4 (the dotted lines). We compared these in the u) = and u) = oo cases. Static solutions can exist to the left side of 
each line. For both solutions, v max becomes large in BD theory. For the self-gravitating monopole, the nontrivial structure is 
larger than black hole's, so the change of Vmax from GR becomes larger than its black hole's. 
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FIG. 5. Schematic diagram of the existence region in v-Th diagram. In the shadowed region, the nontrivial solution exists. 
For a fixed value of ru < f* H , there exists the region (i) (v e xtreme < v < Vmax) where two nontrivial solutions with different 
mass exist. On the other hand, for a fixed value of v < v* , there exists the region (ii) {th,rn < f < rH.max) where two 
nontrivial solutions with different mass exist. The region (i) and (ii) disappears for A > 0.287 and A > 0.696, respectively. 
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FIG. 6. Field configurations of the monopole black hole for A = 0.1, v — 0.1, uj = and the horizon radius fu — 0.2, 0.4, 0.6 
((a)f-/i, w (b) r-fhg (c) r-6 (d) f-(p). We also plot the self-gravitating monopole in the solid line (i.e., fu = 0) for comparison. 
As the horizon radius becomes large, the nontrivial structure becomes small. So the variation of the BD scalar field becomes 
small. 





FIG. 7. Mass-horizon radius diagram of the monopole black hole for v — 0.1, A = 0.1 and u) = —1.4, — 1, 0, oo. ((b) 
is a magnification of (a) around the maximum horizon radius.) The solid lines denote monopole black holes and the dotted 
lines denote RN black holes. Unlike the RN black hole, the monopole black hole does not have an extreme limit but has the 
f H — limit, which corresponds to the self-gravitating monopole solution. The high entropy branch ends up with the maximum 
horizon radius and the maximum mass at point A, where a cusp structure appears. The low entropy branch is connected to 
the RN black hole branch. 
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FIG. 8. Field configurations of the solution in the high entropy branch for th = 0.59, A = 0.1, v — 0.1 and w = —1.4, — 1, 
0, oo ((a) r-h, w (b) r-rh g (c) r-ip ). The YM field and the Higgs field hardly depend on w. 
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FIG. 9. Field configurations of the solution in the low entropy branch for fn — 0.59, A = 0.1, V — 0.1 and u> — —1.4, — 1, 
0, oo ((a) f-h, w (b) f-rh g (c) r-if> ). The YM field and the Higgs field depend on u contrary to the high entropy branch. It is 
because ?h — 0.59 line crosses near the point A for small to, but it crosses near the point B for large to. 
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FIG. 10. lo dependence of (a) the gravitational mass and (b) the scalar mass of the monopole black hole for v — 0.1 and 
A = 0.1. We choose th — 0.59 for which there exist two branches of monopole black holes for fixed to. Though dependence of 
M in both branches can be approximated by a RN black hole, the effect of the nontrivial structure appears clearly for M s . 
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FIG. 11. Mass-inverse temperature diagram for (a) v — 0.1, A = 1 and to — 0, oo, (b) v — 0.1, A = 0.1 and u) = —1.4, — 1, 0, 
oo. When we consider the Hawking evaporation of a RN black hole, the solution experiences the second order phase transition 
for A = 1 when the solution changes to the monopole black hole at the point B. On the contrary, it experiences first order 
phase transition for A — 0.1. 
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FIG. 12. Gravitational mass- horizon radius diagram of the monopole black hole in BD frame for v = 0.1, A = 0.1 and 
u — —1-4, — 1, 0, oo. The solid lines denote the monopole black hole and the dotted lines denote the RN black hole. In this 
diagram, a cusp structure disappears. This may show that Mbd is not appropriate to the control parameter of catastrophe 
theory. 
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